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Section  1 
INTRODUCTION 


It  has  long  been  known  that  the  magnetosphere  contains  many  regions  In  which 
energetic  particles  are  trapped.  The  Van  Allen  belts  are,  of  course,  the 
classical  example  containing  particles  with  energies  up  to  and  beyond  1 Mev. 
There  are,  however,  other  regions  In  the  tall  of  the  magnetosphere  and  In  the 
Inner  magnetosphere  where  additional  relatively  energetic  particle  populations 
are  found.  Several  detailed  explanations  have  been  given  for  the  existence 
of  the  very  high  energy  Van  Allen  radiation.  The  remaining  energetic  particle 
populations,  however,  must  be  explained  In  terms  of  the  Interaction  of  their 
source  region,  the  solar  wind,  with  the  geomagnetic  field.  The  solar  wind  has 
particle  energies  on  the  order  of  1 keV  while  the  energetic  particle  populations 
within  the  magnetosphere  have  energies  which  range  from  several  keV  to  several 
tens  of  keV.  It  Is,  therefore,  necessary  to  explain  the  energization  of  these 
particles  In  terms  of  ’'local"  acceleration  mechanisms.  That  Is,  the  energi- 
zation of  these  particles  must  take  place  within  the  magnetosphere  or  In  the 
blow  shock  and  magnetosheath  regions  which  surround  It.  There  are  only  two 
ways  In  which  charged  particles  can  be  accelerated:  either  by  time  varying 
magnetic  fields  or  the  presence  of  electrostatic  fields. 

Generically  there  are  two  kinds  of  electric  fields:  electrostatic  fields 
result  from  processes  which  cause  charge  separation  over  a region  of  space,  and 


time  varying  magnetic  fields  which  have  associated  induced  electric  fields.  In 
describing  both  types  of  fields,  care  must  be  taken  in  selecting  the  proper 
coordinate  system.  Because  reference  frames  are  so  important  with  electric 
fields,  the  transformation  of  magnetic  and  electric  fields  from  one  reference 
frame  to  another  must  be  considered.  Such  transformations  often  result  in 
the  presence  of  the  electric  field. 

In  the  earth's  ntagnetosphere  several  electric  fields  have  been  postulated  to 
explain  observed  plasma  phenomena.  The  electric  fields  associated  with  the 
rotation  of  the  earth  are  used  to  explain  the  corotation  of  the  plasma  which 
populates  the  plasmasphere.  A "convection"  electric  field  arising  from  some 
unknown  viscous  interaction  between  the  magnetosphere  and  the  solar  wind  has  been 
postulated  to  explain  gross  plasma  motions  in  the  magnetosphere.  More  recently 
electric  fields  in  the  ionosphere-magnetosphere  system  have  been  suggested  to 
explain  the  existence  of  currents  flowing  along  magnetic  field  lines.  Thus, 
the  electric  fields  that  have  been  discussed  in  the  literature  to  date  have 
been  linked  with  specific  observations  in  the  magnetoshere  but  have  not  typically 
addressed  the  general  question  of  particle  acceleration  and  energization. 


It  is  useful  then  to  consider  what  we  know  quantitatively  about  the  existence 
of  electric  fields  in  the  magnetosphere.  It  turns  out  that  little  is  known 
concerning  the  electrostatic  fields  caused  by  charge  separation.  The  subject 
of  induced  electric  fields  has  also  received  little  attention  to  date.  However, 
the  source  of  an  induced  electric  field  must  be  a time  varying  magnetic  field. 
Considerable  work  has  recently  been  done  on  the  nature  of  time  varying  magnetic 


fields  In  the  magnetosphere.  Thus,  the  source  of  at  least  some  Induced 
electric  fields  Is  understood  quantitatively.  It  Is  because  of  this  availa- 
bility of  Information  concerning  the  source  functions  that  we  have  during  the 
past  year  quantitatively  examined  the  Induced  magnetospheric  electric  field. 

From  the  point  of  view  of  particle  energization  Induced  electric  fields  are 
somewhat  more  Interesting  than  electrostatic  fields.  The  electrostatic  field 
Is  conservative.  Thus,  a particle  moving  through  It  upon  return  to  Its 
Initial  position  will  end  up  with  no  net  energy  change  although  It  may  undergo 
energy  changes  along  Its  path.  The  Induced  electric  field,  however.  Is  not 
conservative  and  particles  moving  In  It  can  be  accelerated  and  energized  even 
If  they  move  over  a closed  path.  The  selection  of  a convenient  reference  frame, 
and  Its  Impact  on  the  electric  fields  which  can  be  observed  In  rotating  reference 
systems  are  discussed  elsewhere  In  this  report. 

The  goal  of  this  work  Is  to  ultimately  shed  light  on  the  local  acceleration 
processes  which  must  exist  In  the  magnetosphere  and  that  are  responsible  for 
the  local  energization  of  charged  particles.  During  the  past  year  we  have 
examined  In  quantitative  detail  the  Induced  electric  field  associated  with  the 
dally  wobble  of  the  earth's  magnetic  field  and  the  response  of  the  magnetospheric 
current  systems  to  It.  We  have  found  that  this  source,  which  Is  always  present. 

Is  smaller  than  the  Induced  electric  field  expected  during  substorm  and  storm 
periods.  Although  It  Is  found  to  be  a small  fraction  of  a millivolt  per  meter 
throughout  most  of  the  magnetosphere.  It  still  Is  large  enough  to  be  of 
Interest  In  the  study  of  particle  energization.  The  response  of  the  magnetospheric 
plasma  to  the  presence  of  this  Induced  electric  field  Is  also  discussed.  The 
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response  mechanism  is  very  complicated  and  in  fact  different  mechamisms 
may  be  operative  in  different  regions  of  the  magnetosphere.  Before  discussing 
this  induced  electric  field  and  the  plasma's  response  to  it,  some 
connents  are  made  on  the  quantitative  modeling  of  magnetic  fields  and  in 
particular  the  magnetospheric  magnetic  field  resulting  from  the  daily  wobble 
of  the  geomagnetic  field  as  it  responds  to  the  solar  wind. 


Section  2 


QUANTITATIVE  MODELS  OF  MAGNETOSPHERIC  MAGI^ETIC  FIELDS 


The  authors  of  this  report  have  developed  with  partial  support  from  ONR 
several  models  of  the  magnetospheric  magnetic  field.  These  models  differ 
from  their  predecessors  in  that  they  include  all  three  major  magnetospheric 
current  systems;  magnetopause,  ring,  and  tail.  Because  the  ring  and  tail 
currents  flow  over  a volume  of  the  magnetosphere,  it  was  required  that  a pro- 
cedure be  developed  for  quantitatively  representing  such  distributed  currents. 
Previous  models  had  used  only  point  sources  (dipole  field)  and  currents  flowing 
on  surfaces  (on  the  magnetopause  and  neutral  sheet  regions).  The  details  on 
how  to  represent  these  distributed  currents  have  been  discussed  in  other 
reports  and  published  papers.  It  was  also  required  that  a new  method  for 
representing  the  magnetic  field  be  found.  Previously,  models  had  been  described 
in  terms  of  a spherical  harmonic  series  which  represented  the  scalar  magnetic 
potential.  Derivatives  of  this  potential  yield  the  components  of  the 
magnetic  field.  This  formalism,  however,  cannot  be  used  in  the  prejsence 
of  the  currents  since  the  field  from  the  scalar  potential  is  curl  free.  A 
representation  of  the  field  which  allowed  the  oresence  of  a finite  VxF  was 
required.  Details  of  this  problem  have  also  been  published.  The  resulting 
magnetic  fields  have  been  tested  by  comparing  their  output  with  many  observed 
particle  and  field  features  in  the  magnetosphere.  These  include  comparison 
with  AB  contours  of  Sugiura  and  Porof,  the  correct  calculation  of  observed 
low  latitude  cutoff  boundaries  for  precipitating  solar  cosmic  ray  particles. 


5 


I 


the  high  latitude  boundary  of  the  trapping  region  of  low  energy  particles,  the 
topology  of  the  field  as  compared  with  barium  cloud  releases,  and  detailed 
comparisons  of  model  calculations  with  the  observed  field  measured  at  geosyn- 
chronous orbit.  The  resulting  magnetic  field  topology  is  shown  for  the  noon- 
midnight  meridian  plane  in  Figure  1.  It  is  noted  that  the  field  lines  are  much 
less  dipolar  than  in  previous  models.  This  is  the  result  of  the  inclusion  of  the 
contributions  from  the  distributed  quiet  time  ring  current,  which  because  of 
its  diamagnetic  nature  depresses  the  magnetic  field  in  the  inner  magnetosphere. 

This  was  an  output  from  our  first  model  in  this  series.  It  was  developed  with 
the  restriction  that  the  solar  wind  be  incident  perpendicular  to  the  earth's 
dipole  axis.  In  a more  recent  model  this  restriction  was  relaxed  and  all 

"tilt  angles"  were  permitted.  The  resulting  model  can  be  used  to  represent  the 

daily  and  seasonal  variations  in  the  quiet  time  magnetospheric  magnetic  field 
configuration.  This  model  was  tested  extensively  using  quiet  time  geosynchronous 
magnetometer  data.  It  is  this  model  that  is  of  interest  to  the  present  study 
because  such  a time  varying  magnetic  field  produces  an  electric  field. 

2.1  TIME  VARYING  MAGNETIC  FIELD 

As  pointed  out  by  Hones  and  Bergeson  (1965),  the  existence  of  only  one  time 
varying  magnetic  field  does  not  result  in  the  net  energization  of  charged 
particles.  This  is  because  a reference  frame^an  always  be  chosen  that  moves 
with  the  magnetic  field  variation  such  that  in  that  frame  no  variation  persists. 
Thus,  for  the  case  of  the  wobbling  dipole,  in  a geographic  reference  frame 
(in  which  the  dipole  is  fixed)  3B/9t  will  be  zero.  However,  when  the  response 

of  the  magnetospheric  magnetic  field  is  included,  there  can  be  no  coordinate  system 

in  which  both  the  earth's  dipole  motions  and  the  resultant  changes  in  magnetospheric 
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current  systems  can  be  cancelled  out  such  that  9B/dt  equals  zero.  It  is  thus 
expected  in  the  present  problem  that  a finite  induced  electric  field  will 
result.  It  was  the  goal  of  the  current  periq.d  of  performance  to  describe  this 
induced  electric  field  quantitatively. 

2.2  VECTOR  POTENTIAL  REPRESENTATION 

Maxwell's  equation  7xF  = -3B/3t  suggests  a relation  between  the  time  varying 
magnetic  field  and  the  electric  field  but  does  not  provide  a simple  means  for 
obtaining  the  induced  electric  field  from  the  time  varitions  in  the  magnetic 
field.  It  is  therefore  convenient  mathematically  to  work  with  the  magnetic 
vector  potential. A^.  In  all  of  our  recent  modeling  efforts  we  have  routinely 
determined  If  when  we  calculate  the  magnetic  induction  F.  S'  is  found 
much  the  same  way  as  B by  integrating  over  the  current  system  according  to 
the  formula  S'  = Uq/Att  / JdV/r.  The  use  of  S'  provides  several  advantages. 

It  can  be  used  directly  to  determine  F,  from  F » VxR.  This  is  required  in  cases 
where  V*F  must  be  rigorously  zero.  (Our  usual  model  consists  of  a direct 
representation  for  F and  does  not  require  taking  the  Vxff,  eliminating  additional 
numerical  complications.)  The  time  variation  in  B also  directly  yields  the 
value  of  induced  electric  field.  Thus,  the  field  produced  throughout  the 
magnetosphere  because  of  the  daily  wobble  in  the  earth's  dipole  axis  can  be 
determined  quantitatively.  We  believe  that  this  is  the  first  example  of  an 
electric  field  system  in  the  magnetosphere  being  represented  quantitatively 
with  a firm  understanding  of  its  source. 
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Section  3 

THE  QUANTITATIVE  REPRESENTATION  OF  THE  INDUCED  ELECTRIC  FIELD 


In  the  detenni nation  of  the  magnetic  field  produced  by  the  magnetospheric 
current  systems  it  is  necessary  to  integrate  over  the  currents  using  the 
Biot-Savart  law  to  yield 


f ■ dV 

where  r is  the  vector  from  the  element  of  current  to  the  point  where  is  being 
found.  Because  of  the  complexity  of  this  current  system,  the  Integral  is  performed 
numerically.  Similarly,  and  with  very  little  additional  cost  in  computer  time, 
the  vector  potential  for  this  field  can  also  be  obtained. 


The  induced  component  of  the  magnetospheric  electric  field  resulting  from  this 
dependence  of  the  magnetic  field  explicitly  on  the  tilt  angle  and  thus  on  time 
can  be  obtained  from  the  vector  potential 


(Appendix  A describes  a method  for  calculating  the  vector  potential  for  IT  of 
any  arbitrary  tilt  angle.) 
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I This  procedure  for  obtaining  the  magnetospheric  Induced  electric  field  (and 

the  magnetic  field)  Is  actually  an  approximation.  The  existence  of  a medium 
In  which  these  fields  are  present  Is  considered  only  as  a source  (3)  of  these 
fields;  the  fields  have  traditionally  been  considered  as  virtually  static  and 
all  propoagatlon  effects  Ignored.  In  regions  where  the  plasma  Is  particularly 
tenuous  this  approximation  yields  satisfying  results,  particularly  In  the 

1 case  of  the  magnetic  field.  The  Induced  electric  field,  however,  presents  a 

[ 

i problem.  In  that  In  addition  to  the  vector  potential  contributions,  the  scalar 

potential  component  = -7(J»  must  also  be  considered.  The  total  field  ^ j “ - 
Fj  - -V(j>-dA/3t  Is  the  observable  quantity,  and  both  aspects  of  the  field  must 
be  properly  modeled.  It  Is  certain  that  the  presence  of  an  Ionized  medium 
In  the  magnetoshere  Is  Important  to  3A/dt  but  that  under  certain  conditions, 
discussed  below.  It  can  be  Ignored  to  a first  approximation.  In  the  case  of 
(the  scalar  potential  component)  both  the  local  medium  as  well  as  extended 
regions  contribute  to  It.  Thus,  the  total  electric  field  at  a given  location 
will  not  only  Include  the  contribution  from  the  local  time  dependence  of  the 
magnetic  field,  but  also  contributions  from  other  regions  of  the  magnetosphere 
transmitted  to  this  location  by  the  medium,  to  the  extent  that  It  Is  conductive. 

3.1  THE  SCALAR  POTENTIAL  ELECTRIC  FIELD  IN  PRESENCE  OF  AN  IONIZED  MEDIUM 
A procedure  developed  by  Hones  and  Bergeson  (1965)  \ias  applied  to  evaluate  the 
scalar  potential  component  of  the  electric  field.  They  assumed  that  In  the  region 
of  Interest,  charged  pa<^t1cle  densities  are  high  enough  to  assure  sufficient  con- 
ductivities along  the  magnetic  field  lines  such  that  Tj*!  ■ 0. 

7i,-_  (W  8A„ 
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where  d4/ds  Is  the  scalar  potential  gradient  along  the  field  line,  and  dA„/dt  Is 
the  component  of  Tj  parallel  to  the  magnetic  field  line.  This  equation  can  be 
Integrated  along  thi>  field  line  to  yield  the  scalar  potential.  Thus, 


* ♦b  ■ ^ * 

where  the  line  Integral  is  along  the  field  line,  and  Is  the  potential  at  some 
appropriate  boundary. 

Because  this  approach  greatly  simplifies  the  problem.  It  Is  ten^tlng  to  use  It,  along 
with  models  of  3A/3t,  to  quantitatively  determine  possible  sources  of  the  electric  field 
observed  In  the  magnetosphere.  The  presence  of  the  ionized  medium,  which  makes 
this  approach  possible,  however,  also  complicates  the  problem.  Propagation 
effects  on  iJ/it  must  be  considered,  as  well  as  the  shielding  of  the  contributions 
which  depend  strongly  on  the  selection  of  the  appropriate  boundary  condition. 

The  effects  of  the  medium  on  the  field  line  Integral  of  3A/3t  are  summarized  In 
Table  1.  The  simple  case  of  time  Independent  B Is  included  here  only  for  com- 
pleteness, since  It  does  not  represent  a very  Interesting  example.  It  Is  consis- 
tent, however,  with  the  assumption  of  magnetic  field  lines  being  equi potentials  In 
an  Ionized  medium.  The  second  case  represents  the  Hones  and  Bergeson  procedure 
(discussed  above)  which  is  appropriate  In  regions  where  the  conductivity  along 
field  lines  Is  high  enough  so  that  - 0,  but  low  enough  to  allow  nonlocal 
contributions  to  the  total  electric  field.  In  the  extreme  case  of  high  conduc- 
^ tivltles  (Case  3)  along  field  lines,  only  locally  Induced  electric  fields  need 
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CONDITIONS  FOR  SELECTIONS  APPROPRIATE  VALUES  OF  THE  INDUCED  ELECTRIC  FIELD 


be  considered;  external  sources  can  contribute  only  to  the  extent  that  the  fields 
originating  from  them  can  penetrate  the  medium  without  energy  dissipation. 


Case  4,  In  Table  1.  represents  the  conditions  present  In  areas  whre  the  conductivity 
along  field  lines  Is  so  low  that  T-F  ^ 0.  The  procedures  discussed  here  for  the 
solution  of  I are  then  not  appropriate  and  the  only  available  approach  Is  the 
solution  to  the  wave  equation.  This,  however.  Is  the  vacuum  case  and  probably 
does  not  apply  anywhere  In  the  magnetosphere. 

3.2  BOUNDARY  CONDITIONS 

The  Hones  and  Bergeson  procedure  can  be  applied  to  large  regions  of  the  magnetosphere 
to  obtain  first  order  estimates  of  the  total  electric  field.  Contributions  from  those 
regions  of  the  magnetosphere  where  this  procedure  Is  Inappropriate,  can  be  approxi- 
mated by  the  selection  of  other  appropriate  boundary  conditions.  Models  of  Ion 
concentrations  throughout  the  magnetosphere  are  required  for  this  purpose.  When 
these  are  not  available,  limited  areas  of  the  magnetosphere  can  be  Investigated 
using  available  experimental  data. 

Other  regions  of  the  magnetosphere  can  be  Investigated,  by  utilizing  the  fact 
that  the  magnetic  field  lines  terminate  In  the  Ionosphere.  Neglecting  the  presence 
of  high  conductivity  regions  outside  the  Ionosphere,  the  Ionosphere  can  be  considered 
f as  the  boundary.  The  value  of  the  scalar  potential  In  the  Ionosphere  Is  still 

1 subject  to  speculation.  One  candidate  based  on  theoretical  consideration  Is  the 

potential  of  a unlformally  magnetized  rotating  sphere  which  was  used  by  Hones  and 
Bergeson  (1965)  In  their  calculations.  This  potential  Is  given  by 
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where  V » the  potential  at  point  r, 

U » w X r , 

» W X r/R^, 

w « angular  velocity 
y a magnetic  moment 
> radius  of  sphere. 

Another  candidate  boundary  condition  would  be  to  simply  consider  the  upper 
Ionosphere  as  an  equipotentlal  region.  This  assumption,  as  well  as  the  uni- 
formal  1y  magnetized  sphere  model  were  employed  In  the  present  work  to  evaluate 
the  sensitivity  of  the  electric  field  model  to  both  the  choice  of  the  boundary 
condition  and  the  nonlocal  contributions  to  the  Induced  electric  field.  Some 
results  of  this  Investigation  are  Illustrated  In  Figures  2 and  3.  In  both 
figures,  equipotentlal  contours  are  plotted  In  the  geographic  equatorial  plane. 
The  sun  Is  to  the  left  (negative  x-directlon)  In  both  plots.  The  calculations 
were  done  In  an  Inertial  coordinate  system  In  which  the  earth  is  rotating,  the 
Induced  electric  field  (i.e.,  dA„/3t)  Is  due  to  both  the  dipole  field  and  the 
diurnal  variations  In  the  current  systems.  In  both  figures  four  calculation 
times  are  Illustrated:  thus,  for  example.  In  Figure  2a  the  calculation  time 
Is  UT  (universal  time)  - 0.0  hours,  summer  solstice,  when  the  Greenwich  meridian 
(positive  x-directlon)  Is  at  midnight.  In  Figure  2b  the  calculation  Is  at 
UT  - 6.0  hours  when  the  Greenwich  meridian  Is  at  dawn  (I.e.,  In  the  positive 
y-direction).  Similarly,  Figures  2c  and  2d  are  at  UT  - 12.0  and  18.0  hours 
when  the  Greenwich  meridian  Is  at  noon  (the  negative  x-directlon)  and  dusk 
(the  negative  y-directlon),  respectively. 
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In  Figure  2 a uniform  (constant)  potential  surface  at  r = 1 was  used  as  the 
boundary  condition;  while  in  Figure  3 the  rotating,  uniformly  magnetized  shphere 
was  used  at  the  boundary  of  r = 1 R^.  The  sensitivity  to  boundary  conditions 
is  clearly  evident.  The  asymmetry  in  both  sets  of  figures  is  due  to  the 
induced  electric  field  resulting  from  time  variation  of  the  magnetosheric 
current  systems.  This  is  particularly  noticeable  in  Figure  3 where  without 
the  current  system  the  contours  would  be  circular. 

In  Figure  3 the  change  in  potential,  ((>,  between  the  ionosphere  (where  it  is 
assumed  to  be  zero)  and  the  equatorial  plane  is  illustrated.  It  is  evident 
from  this  figure  that  if  the  ionospheric  electric  field  (and  ()))  is  variable, 
the  total  electric  field  in  the  magnetoshere  will  also  vary  substantially  with 
time.  Thus,  through  this  mechanism,  the  ionosphere  may  exert  some  influence 
on  the  dynamics  of  magnetos pheric  fields  and  low  energy  plasma.  Also,  Figure  2 
indicates  the  change  in  <})  between  the  ionosphere  and  the  equatorial  plane  when 
the  ionospheric  value  of  (|)  is  taken  to  be  zero.  Thus,  magnetic  field  lines  are 
not  equipotentials  even  when  only  the  "small"  induced  fields  like  the  wobbling 
dipole  are  present. 

The  various  components  of  the  total  electric  field  which  would  be  "felt"  by  a 
charged  particle  rotating  with  the  earth  in  geosynchronous  orbit  are  illus- 
trated in  Figures  4 and  5.  The  uni  formally  magnetized,  rotating  sphere  boundary 
condition  was  used  for  these  calculations.  Figure  4 displays  the  field  as  a 

V 

function  of  longitude;  Figure  5,  the  electric  field  as  a function  of  time. 
Contours  of  constant  the  equatorial  plane  are  shown  in  Figures  6 
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Figure 


Total  Electric  Field  Experienced  by  a Charged  Particle  Rotating 
with  Earth  at  Synchronous  Orbit 


through  8 for  the  rotating  magnetized  sphere  boundary  condition.  Figure  6 
illustrates  the  total  electric  field  in  an  inertial  frame  at  four  times 
(UT  = 0,  6,  12,  and  18  hours,  respectively).  The  contours  in  Figure  6,  like 
those  in  Figures  2 and  3,  are  plotted  in  inertial  coordinates  which  coincide 
with  the  geographic  coordinates  at  UT  = 0 hours;  the  sun  is  to  the  left  (in 
the  negative  x-direction) . In  Figure  7 the  calculations  illustrated  were  done 
in  the  geographic  (i.e.,  rotating)  coordinate  system;  however,  the  plots  are 
displayed  such  that  the  location  of  the  sun  is  the  same  as  in  the  inertial 
coordinate  plots.  The  location  of  Greenwich  is  indicated  by  the  arrow  in 
Figure  7.  In  Figure  8 the  constant  magnetitude  contours  for  l^total  ^ 
rotating  system  are  shown  for  the  alternate  boundary  condition  of  <|)g(^Q)  = 0 
at  r„  = 1 Rn. 

The  model  of  the  total  electric  field  discussed  above  contains  contributions 
from  t.je  time  varying  magnetic  field  as  well  as  the  electrostatic  contributions 
which  would  be  contained  in  the  boundary  condition  on  the  scalar  potential.  The 
sensitivity  of  this  electric  field  model  to  the  boundary  conditions  has  been 
discussed  above.  The  induced  electric  field  contributes  to  the  model  in  two 
ways:  in  the  direct  contributions  of  the  induced  electric  field, 

= -dl/dt, 

and  in  the  form  of  the  line  integral  contribution  to  the  scalar  potential 
described  above.  The  direct  contribution  of  the  induced  electric  field  to  the 
total  electric  field  illustrated  in  Figure  6 is  given  in  Figure  9.  All 
comments  about  the  coordinate  system  made  with  respect  to  Figure  6,  also 
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Figure  8.  I^totall  ® rotating  frame  for  a uniform  boundary  condition. 
(Units  are  volts/meter  and  the  major  divisions  are  5 R^.) 
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Figure  9.  |Fj|  in  the  Inertial  coordinate  system.  (Units  are  volts/meter, 
the  major  divisions  are  5 R^.) 
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apply  to  Figure  9.  The  range  of  values  for  |Ej|  Illustrated  are  1.0  to 
0.01  mV/m,  although  only  portions  of  the  1 mV/m  contour  can  be  seen  In  the 
figure  (the  1.0  mV/m  contour  Is  nearest  the  1 R^,  earth  outline).  In  regions 
where  the  Induced  field  Is  a maximum.  Its  direct  contribution  Is  seen  to  be 
about  14  percent  of  the  total  field.  Since  the  Induced  electric  field  Is 
presented  In  Inertial  coordinates,  the  dipole  contributions  are  seen  to  domi- 
nate for  small  values  of  R.  Distortions  from  the  dipole  pattern  are  due  to 
the  diurnal  variations  In  the  external  current  system  (although  these  are  also 
due  to  the  motion  of  the  dipole);  these  are  seen  to  become  significant  at 
R > 4 R^  In  regions  of  maximum  |Ej|,  but  are  Important  over  most  values,  R, 

In  the  region  of  minimum  |Ej|  (I.e.,  In  the  cusps). 

3.3  CHARGED  PARTICLE  ACCELERATION 

To  Indicate  how  fields  such  as  those  Illustrated  In  Figures  4 and  5 would 
contribute  to  charged  particle  acceleration,  the  time  rate  of  change  of  a 
proton's  kinetic  energy  was  calculated  using  the  formula  from  Northrop  (1963) 

-31-  “ q E-V  + M ^ , 

where  <w>  Is  the  kinetic  energy  averaged  over  a gyration  period,  M Is  the  magnetic 
moment  of  the  proton,  q Its  charge,  and  V the  velocity  of  Its  guiding  center. 

The  results  of  this  calculation  are  shown  In  Figures  10  and  11.  Figure  9 shows 
the  energy  gained  by  the  guiding  center  of  motion  as  a function  of  the  guiding 
center  kinetic  energy,  and  similarly.  Figure  10  shows  the  energy  change  In  the 
motion  about  the  guiding  center. 

C 


25 


Figure  Time  Rate  of  Change  of  Proton  Kinetic  Energy,  Guiding  Center  Motion  Contribution 

(UT  - 16.6  Mrs.,  Day  • 172,  Equatorial  Plane,  Longitude  ■ 0*. 


The  energy  gained  by  the  guiding  center  motion  (Figure  10)  corresponds  to  the 
Incremental  change  In  potential  as  a result  of  this  motion.  These  would 
thus  correspond  to  circular  motions  of  the  appropriate  radius  through  potential 
differences  such  as  those  Illustrated  In  Figure  3 but  also  Including  contribu- 
tions from  3B/3t  since  the  total  field  Is  used  In  the  above  expressions.  The 
energy  change  In  the  motion  about  the  guiding  center  (Figure  11)  Is  seen  to 
depend  only  on  3F/3t  and  represents  the  particles  betatron  acceleration. 
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Section  4 


SUMMARY  AND  CONCLUSIONS 

It  Is  an  observed  fact  that  many  regions  of  the  magnetosphere  contain  charged 
particles  with  energies  of  several  tens  of  kilovolts.  It  Is  also  known  that 
the  source  of  these  particles  Is  either  the  earth's  upper  atmosphere  or  the 
solar  wind,  both  of  which  have  particle  enegles  of  one  kilovolt  (the  solar 
wind),  or  only  a few  electron  volts  (the  high  latitude  upper  atmosphere). 

Thus,  It  Is  necessary  to  explain  the  energization  of  charged  particles  occurring 
locally,  that  Is,  within  the  magnetosphere.  Only  electric  and  time  varying 
magnetic  fields  can  produce  this  energization.  It  Is  therefore  Important  to 
understand  and  quantitatively  describe  such  fields  In  the  earth's  magnetosphere. 

Of  the  various  sources  of  electric  fields,  the  time  varying  magnetic  field  has 
been  documented  most  thoroughly.  In  our  work  on  magnetic  field  modeling  we 
have  routinely  computed  the  vector  magnetic  potential  simultaneously  with  the 
calculation  of  magnetic  Induction.  A knowledge  of  the  time  history  of  the 
magnetic  potential  readily  yields  the  associated  Induced  electric  field.  Thus, 
the  Induced  electric  field  from  any  time  varying  magnetic  field  can  be  determined 
If  the  magnetic  field  and  Its  associated  magnetic  vector  potential  are  known 
I quantitatively.  In  the  present  work  we  have  dealt  primarily  with  the  electric 

field  Induced  by  the  dally  wobble  of  the  earth's  dipole  field  and  the  associated 
response  of  the  magnetospheric  current  systems.  Such  a field  Is  Important  ' 

because  It  Is  always  present.  The  results  of  this  work  suggest  quantitatively 
the  presence  of  an  electric  field  that  should  exhibit  some  dependence  on  both  local 
and  universal  time  (because  the  form  of  the  field  Is  dependent  on  the  location  of 
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the  dipole  axis  and  magnetospheric  magnetic  field  topology).  It  was  found 


quantitatively  that  the  raw  Induced  electric  field  is  on  the  order  of  a small 
fraction  of  a millivolt/meter  throughout  most  of  the  magnetosphere.  This 
suggests  that  the  Induced  electric  field  associated  with  the  dally  wobble  of 
the  earth's  magnetic  field  Is  large  enough  to  be  Important  to  the  dynamics  of 
charged  particles  In  the  magnetosphere. 

What  Is  actually  needed  for  the  study  of  particle  field  Interactions  Is  the 
total  electric  field  produced  as  the  magnetospheric  plasma  responds  to  this 
induced  electric  field.  As  might  be  expected,  the  response  of  the  plasma  to 
this  Induced  electric  field,  or  any  other  Induced  electric  field,  depends  on 
both  the  intensity  of  the  magnetic  field  and  various  plasma  parameters.  Our 
analyses  suggest  that  it  may  be  necessary  to  treat  the  response  differently 
in  different  regions  of  the  magnetosphere.  Generally,  however,  it  may  be 
argued  that  the  plasma  responds  to  the  Induced  electric  field  by  tending  to 
short  out  the  field  in  the  direction  of  the  magnetic  field.  This  is  because 
of  the  rather  high  conductivity  that  the  charged  particles  have  along  magnetic 
field  lines.  The  cancelling  out  of  the  induced  electric  field  along  the  direction 
of  the  magnetic  field  is  actually  carried  out  by  the  setting  up  and  maintenance 
of  a very  weak  charge  separation  system.  This  electric  field  can  be  described 
as  a divergence  of  a scalar  potential.  The  problem  of  determining  the  total 
electric  field  then  becomes  one  of  finding  the  appropriate  value  of  this 
scalar  potential.  This  necessitates  the  determination  of  the  boundary  condition 
on  the  scalar  potential.  The  study  of  the  appropriate  boundary  conditions 
on  the  scalar  potential  Is  probably  worthy  of  Its  own  separate  effort.  Our 
work  on  this  subject  to  date  suggests  two  Important  findings;  first  It  Is  not 
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appropriate  to  treat  magnetic  field  lines  as  electrical  equi potentials,  thus 


the  mapping  of  the  Ionospheric  electric  field  to  magnetos pher 1c  regions  must 
Involve  a more  detailed  process.  The  change  In  potential  along  a magnetic 
field  line  In  the  presence  of  the  electric  field  Induced  by  the  wobbling  dipol 
Is  as  large  as  13  kilovolts,  between  the  earth's  surface  and  geosynchronous 
altitudes  (along  a field  line).  The  other  finding  Is  based  on  the  dependence 
of  the  electric  field  on  the  value  of  the  scalar  potential  In  the  Ionosphere. 
This  suggests  strongly  that  as  the  Ionospheric  electric  field  (and  Its 
associated  scalar  potential)  vary  appreciably  as  It  Is  "driven"  by  these 
changes  In  the  upper  atmosphere  and  Ionosphere.  This  suggests  that  some 
magnetospheric  dynamic  processes  may  be  the  result  of  changes  occurring  In 
the  Ionosphere. 


A detailed  examination  of  the  total  electric  field  produced  by  the  wobbling 
dipole  suggests  that  Its  impact  on  particle  energization  and  acceleration  may 
be  an  Important  feature  In  the  vicinity  of  geosynchronous  orbit  during  quiet 
magnetic  conditions.  It  might  be  expected  that  this  field,  since  It  rotates 
somewhat  In  phase  with  the  dipole,  would  have  the  largest  Impact  on  particles 
moving  around  the  earth  with  the  same  angular  velocity.  It  turns  out  that 
one  kilovolt  protons  (the  solar  wind  - the  raw  material  for  the  magnetospheric 
plasma  Is  comprised  of  approximately  one  kilovolt  protons)  take  about 
one  day  to  move  around  the  earth  at  geosynchronous  altitude.  Thus,  the 
Induced  electric  field  may  be  an  Important  mechanism  for  energizing  some 
of  the  low  energy  particles  as  they  enter  from  outside  of  geosynchronous 
orbit  and  move  radially  Inward  through  the  outer  Van  Allen  region. 
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The  actual  energization  process  has  not  yet  been  studied.  However, 

It  Is  known  that  this  particular  electric  field  Is  not  conservative  and  can 
permanently  energize  charged  particles.  Also,  It  would  be  expected  because  of 
phase  differences  between  the  electric  field  and  the  charged  particle  motions 
at  different  locations  In  geosynchronous  longitude,  that  only  a fraction  of  the 
particles  would  actually  be  energized  and  others  may  lose  energy.  However,  It 
Is  well  known  that  even  In  the  outer  Van  Allen  zone  the  energized  particles 
; comprised  only  a small  fraction  of  the  total  plasma  density.  The  effect  of 

this  mechanism  on  charged  particles  then  might  be  summarized  by  saying,  "You 
can  accelerate  some  of  the  particles  some  of  the  time  - and  that  may  be  enough." 

i 

i 
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It  Is  believed  that  this  study  constitutes  the  first  quantitative  examination  of 
an  electric  field  in  the  magentosphere  where  its  source  is  understood  quantitatively. 
The  particular  mechanism  examined,  that  produced  by  the  dally  wobble  of  the  dipole 
magnetic  field  and  the  response  of  the  magnetosphere  to  It,  produces  an  Induced 
electric  field  that  Is  present  all  the  time,  and  large  enought  to  be  of  Interest 
to  the  study  of  charged  particle  dynamics.  It  is  hoped  that  as  new  observational 
electric  field  data  become  available,  that  some  effort  be  made  to  search  for 
local  time  effects  in  the  data.  We  would  hope.  In  the  near  future,  to  do  two 
things:  first,  to  study  the  boundary  conditions  on  the  scalar  potential  for 
the  electric  field  In  more  detail  as  It  is  determined  by  Ionospheric  electric 
fields;  and,  secondly,  to  quantitatively  describe  an  electric  field  Induced  by 
time  varying  magnetic  fields  associated  with  the  magnetospheric  substorm  and 
magnetic  storm  processes.  Finally,  the  response  of  particles  and  plasma  to 
such  combined  electric  and  magnetic  fields,  which  have  been  described  quantitatively 
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in  a time  varying  setting,  should  be  studied  in  detail.  Hopefully  such  self 
consistent  calculations  will  reveal  a more  detailed  understanding  of  several 
dynamic  magnetospheric  processes. 
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Appendix  A 

ELECTRIC  FIELD  CALCULATIONS 

A computer  program  was  developed  to  perform  the  electric  field  calculations 
described  in  this  report.  The  key  to  these  calculations  are  two  subroutines 
(BXYZMU  and  AXYZ)  which  contain  a model  (Olson-Pfitzer)  of  the  earth's 
external  magnetic  field,  B^,  and  the  associated  vector  potential, 

(respectively).  (The  Olson-Pfitzer  model  is  described  in  Reference  Al,  and 
the  model  for  the  vector  potential  in  Appendix  B of  this  report.) 

All  electric  field  calculations  are  performed  in  an  earth  centered,  inertial 
coordinate  system.  Motion  of  the  earth  around  the  sun  is  important  only  to 
the  extent  if  affects  the  external  magnetic  field  (i.e.,  the  tilt  angle 
dependence)  otherwise  it  is  insignificant  and  is  neglected  in  these  calculations. 
Because  of  this  choice  of  coordinate  system,  the  total  magnetic  induction 
vector.  By,  and  vector  potential.  Ay,  must  include  contributions  from  the 
earth's  dipole  field.  Thus,  to  the  B^  and  A^,  obtained  from  subroutine 
BXYZMU  and  AXYZ,  are  added 


A A 


where  p is  the  dipole  moment 
^ r 

and  n * - . 
r 
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The  Induced  electric  field  Fj  = -3S/3t  is  obtained  by  two  calls  to  AXYZ:  one 
at  time  t and  the  other  for  time  t+At.  This  gives  3Ag/3t  which  must  be  combined 
with 


where  u is  the  dipoles  angular  velocity. 

The  total  electric  field  which  is  generated  by  the  temporal  variations  of  the 
magnetic  field  is  expressed  by 


The  scalar  potential,  , is  obtained  from  the  integral 

♦(-■I  ■ 

as  described  in  this  report.  This  integral  is  evaluated  numerically  (using  the 
S.  Gill  technique  [Gill,  1951]).  To  obtain  the  gradient  of  4>,  <p(r)  must  be 
evaluated  at  four  points:  these  are  selected  to  be  r,  r + §ids,  r + 6§2* 
r + Afij,  where  ds  is  the  integration  step,  and  6 an  appropriate  incremental 
step  size.  The  unit  vector  is  tangent  to  Bj  at  point  r,  62  is  a unit  vector 
parallel  to  the  principal  normal  (obtained  from  d6^/ds),  and  §2  “ x 62- 
Thus, 
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(|»(r+ds  §,)  - (|)(r) 

= 3i 


®1  " 


((•(r+6e2)  - 


4»(r+6e3)  - (|)(r) 


The  direction  of  integration  (i.e.,  of  6^)  is  chosen  so  that  the  integral 
terminates  at  r - 1 Rr.  where  a boundary  value,  (|)o  « is  evaluated. 

0 C 0 0 

Two  options  are  currently  available  in  this  program 


♦(^o)  = 0, 

a uniform  boundary  condition,  and 


ioixr^)  (pxr^j) 


the  uniformally  magnetized  rotating  sphere  boundary  condition. 


The  above  procedure  yields  the  total  electric  field. 


in  the  inertial  coordinate  system.  A particle  rotating  with  the  earth  will 
experience  a net  electromagnetic  force  F ■ q(Ej  + v x Bj),  where  v ■ ojxr,  and 
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q is  the  charge  on  the  particle.  This  force  transforms  as  a simple  vector 
from  the  inertial  coordinate  system  to  a coordinate  system  rotating  with  the 
earth.  Thus,  the  electric  field  observed  by  this  rotating  particle  is  given  by 


^T,rot 


where  both  vector  quantities  are  expressed  in  the  rotating  coordinate  system. 


The  two  magnetic  field/vector  potential  routines,  BSYZNU  and  AXYZ,  are  formulated 
in  yet  another  coordinate  system:  The  solar-magnetic  coordinate  system  in 
which  the  Z-axis  is  antiparallel  to  u,  the  earth's  dipole  moment,  and  the 
earth-sun  direction  is  in  the  XZ-p1ane.  This  coordinate  system  was  selected 
for  these  models  (i.e.,  and  because  it  simplifies  their  formulation, 
but  it  does  create  the  necessity  of  numerous  coordinate  transformations  during 
the  calculations.  For  convenience,  a special  subroutine  (COORD)  was  written 
to  perfrom  these  transformations.  Another  subroutine,  ANGLE,  calculates  the 
position  of  the  sun  in  the  celestrial  sphere  and  determines  the  tilt  angle  (u, 
used  in  AXYZ  and  BXYZMU)  as  well  as  the  rotation  sines  and  cosines  for 
transforming  from  geomagnetic  to  solar  magnetic  coordinates. 


Appendix  B 

VECTOR  POTENTIAL  SUBROUTINE  AXYZ 


A listing  of  the  subroutine  to  calculate  the  vector  potential  Is  Included  In 
this  appendix.  This  subroutine  Is  not  In  a final  ANSI  standard  form.  It  Is  the 
Initial  research  code  generated  by  the  least  squares  fitting  routine.  Thus, 
the  FORTRAN  Is  not  necessarily  standard,  efficient  or  easily  transferable  to 
to  other  machines.  The  code  was  written  to  run  on  a CDC  cyber  74  using  the  FTN 
compiler. 

The  routine  has  proven  very  useful  In  calculating  the  Induction  electric 
field.  Since  the  least  squares  routine  generates  the  FORTRAN  directly,  various 
fitting  schemes  were  easily  evaluated. 

The  routine  returns  the  vector  potential  In  units  of  gauss  - R^  for  a given 
location.  The  coordinate  system  used  Is  the  solar  magnetic.  The  calling  sequence 
for  the  routine  Is: 

Input 

XX  A three  dimensional  array  of  the  position  where  the  field  Is 

to  be  calculated. 

XX(1)  Is  In  the  solar  direction  perpendicular  to  X(3). 

XX(3)  Is  parallel  to  the  north  pointing  dipole  axis. 

XX(2)  make  a right  handed  coordinate  system. 

TILT  The  tilt  angle  of  the  dipole  axis  Is  to  the  sun  earth  line. 
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I 


of  the  vector  potential  In  units  of  gauss  - 
ANA6  The  magnitude  of  the  vector  potential. 

A separate  entry  point  to  the  routine,  DAXYZ,  calculates  dA/d(x,  the  derivative 
of  the  vector  potential  with  respect  to  the  t11t  angle  (caution  - entry  points 
are  highly  machine  dependent,  this  feature  for  example  would  not  work  on  an 
IBM  machine).  Note  that  all  coefficients  In  the  listing  were  generated  by 
machine  so  that  no  keypunch  errors  are  present. 
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